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Abstract 

The coexistence of antiferromagnetism with superconductivity is studied theoretically within 
the t-J model with the Zeeman term included. The strong electron correlations are accounted for 
by means of the extended Gutzwiller projection method within a statistically-consistent approach 
proposed recently. The phase diagram on the band filling - magnetic field plane is shown, and sub- 
sequently the system properties are analyzed for the fixed band filling n = 0.97. In this regime, the 
results reflect principal qualitative features observed recently in selected heavy fermion systems. 
Namely, (i) with the increasing magnetic field the system evolves from coexisting antiferromagnetic- 
superconducting phase, through antiferromagnetic phase, towards polarized paramagnetic state, 
and (ii) the onset of superconducting order suppresses partly the staggered moment. The supercon- 
ducting gap has both the spin-singlet and the staggered-triplet components, a direct consequence 
of a coexistence of the superconducting state with antiferromagnetism. 

PACS numbers: 71.27.+a, 74.20.-z, 74.70.Tx, 74.25.Ha 
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I. INTRODUCTION 



The interplay of antiferromegnetism (AF) with superconductivity (SC) is one of the 
important topics in condensed-matter physics,- as better understanding of this subject 
would improve our knowledge of a number of systems such as high-Tc,- heavy-fermion,- and 
organic^ superconductors. In all those systems, superconductivity appears in the vicinity of 
magnetic phases (mostly antiferromagnetic, but also ferromagnetio^"^) . Moreover, magnetic 
interactions or fluctuations are very frequently considered to be the pairing mechanism 
in unconventional superconductors.-"- Typically, antiferromagnetism and superconductiv- 
ity are competing quantum phenomena because of the competition between the Meissner- 
supercurrent screening and the internal-field generation by magnetic ordering. This antago- 
nism can be overcome by a spatial separation of the AF and the SC phases or by subdivision 
of the / electrons into more localized (resulting in AF) and more itinerant (participating 
in SC) parts. However, especially interesting is the situation, when the same electrons are 
involved in both phenomena, as is the case for some heavy-fermion systems. There, SC and 
AF can coexist easily, when the periodicity of magnetic structure Xap{= 2a) is much smaller 
than the coherence length ^ for the Cooper pair. In other words, when a, the staggered 
exchange field averages out to zero within the coherence volume. In this respect the Ce- 
based '115' heavy-fermion compounds - the family of CeMIns (with M = Co, Rh, Ir)^— is 
the most promising, as both antiferromagnetism and superconductivity are believed to arise 
from 4/ electrons, where even the interplay of the two orders can be studied by tuning the 
system with pressure, magnetic field, or doping. 

Also, recently a competitive coexistence of AF and SC has been reported in both 
CeRhln^^""— and CeCo(Ini_a,Cda;)5.— In the latter system, a mutual influence of AF 
and SC has been observed. Namely, it turns out that the onset of SC order with lower- 
ing temperature prevents any further increase of the antiferromagnetic magnetization.— A 
similar type of coexistence has also been observed in CeRhSis.— 

Generally, in the heavy-fermion systems strong correlations among electrons are the rea- 
son for an emergence of new and nontrivial physics. Those nontrivial features should be 
properly accounted for when modeling those systems. In the present paper, an investigation 
of the coexistence of AF with SC in an applied magnetic field is presented. To account for 
strong electron correlations, the Gutzwiller-projected t-J model is used with the Zeeman 
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term included. The extended Gutzwiller scheme proposed recently^ is utihzed for calcula- 
tion of statistical averages of the relevant operators. Our model, although at first sight seems 
too simplified to be related to heavy fermion systems, it nonetheless reflects qualitatively 
principal features observed recently in selected heavy fermion systems. 

It is commonly believed that the minimal model for investigation of heavy-fermion sys- 
tems should be the two-band Periodic Anderson Model (PAM) (see e.g. Ref. l2ll ) or the 
Kondo lattice model.— On the other hand, the one-band calculations have already proved 
fruitful in the analysis of AF and SC coexistence in CeRhlus,— as well as in investigations of 
the high- field low-temperature unconventional superconducting phase of CeCoIns.— The 



narrow-band limit of PAM has been discussed theoretically also elsewhere (see Refs. 



26 



and 



271 . Appendix A). Generally, it appears when only a single hybridized band is involved and in 



the heavy-fermion limit (i.e. when /-level occupancy nj = 1 — 6, with 6 <^ 1). Simply put, 
the t- J- type model reflects the physics of those hybridized and strongly correlated systems 
in the narrow /-band limit. 

The paper is organized as follows. In Sec. |TT] we present the general theoretical formu- 
lation. In Sec. mil we show the numerical results, and finally in Sec. IIVI our findings are 
summarized and an outlook is provided. 



II. MODEL, ORDER PARAMETERS, AND CONSTRAINTS 

We start from the t-J model with the Zeeman term included, as represented by the 
Hamiltonian2^ 



n 



tj = P(^^tijcl^Cj„ + J J^SiSj - h^ahi^^P, (1) 

ija (ij) ia 

where (ij) denotes the summation over bonds, and a = ±1 is the spin z-component. Since 
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for 



its derivation,—!^ the t-J model represents an active field of research (see e.g. Ref. 
a recent analysis of the one-dimensional situation). The t-J model captures the essential 
ingredients of physics of the high-Tc superconductors. The advantage of using this model 
is that both AF and SC come from a microscopic parameter - antiferromagnetic exchange 
J and therefore there are no phenomenological terms in the Hamiltonian (as opposed to 
some earlier studies of AF and SC coexistence). We neglect the orbital effects, as the 
Maki parameter— in the systems of our interest here is high.— The Gutzwiller projector 



P = — fii^fii]^ eliminates double occupancies in real space. In the following we will use 
the more general correlator 

ffc = n.A^^/^Aj-^/^(l-n.,n4), (2) 

where Ajo- are the so-called fugacity factors. Also, this correlator connects the correlated |\E') 
and uncorrelated |\E'o) wave functions,— via 

\^) = Pc\^,). (3) 

This allows to express average of any operator O in the correlated state as 

(6) = (*|6|^) = (4) 

{PcPc), 

where (...)o = (^o|---|^o)- With the above equation one can in principle calculate average 
value of Hamiltonian ([1]) , namely 

W ^ {Hu) = $^t..(cU.) + jJ2(iS!S^) + {S:SJ + SfS])) -hJ2^{n...), (5) 

ija (ij) ia 

but this is a nontrivial task, as after applying the Wick theorem too many terms appear 



(see Ref. |20|, e.g. Eq. (8) and the discussion afterwards), and one has to resort to making 
approximations at this point. There are a few ways to perform this operation, and this is 
still an active field of research, so one can expect new calculation schemes to appear—. Here, 
we use the scheme proposed recently by Fukushima^^i^^ in the local-constraint version, which 
assumes that the average number of particles at any site and with any spin is unchanged by 
the projection, 

{flia) = {nia)o- (6) 

This formalism is known to reproduce the Variational Monte Carlo results better than the 
conventional Gutzwiller approximation (at least, for the projected uniform nonmagnetic d- 



wave BCS superconductor - see Figs. 3 and 4 of Ref. |20| ). The local-constraint version of 
the formalism is quite general in the sense that it is capable of accounting for antiferromag- 
netism, superconductivity, and the ferromagnetic polarization. The explicit expressions for 



all averages appearing in Eq. ([5]) are given in Ref. |20|. To express them in terms of mean- 
fields of our interest, we need to assume what is the character of the uncorrelated wave 
function |\&o)- Since our goal is the description of coexistence of AF and SC, we assume the 



corresponding mean- fields as nonzero at the level of |\l/o) as in the following. We start with 
the particle number in the form 



Hi 



{nia)o = ^(n + a rupM + tuaf e*^""') , (7) 



where n is the band filling (assumed as constant), mpM is the ferromagnetic (longitudinal) 
spin-polarization component, and niAF is the antiferromagnetic (staggered) spin polariza- 
tion. The factor e*^""' (with Q = (7r,7r)) is responsible for the sign reversal of the staggered 
magnetic moment when exchanging the two sublattices A and B.— We also assume the 
superconducting order parameter can be decomposed into two components 

{r^j A^, for i E A — sublattice^ 
(8) 
TijAB, for i G 5 — sublattice, 

where Tij ensures the ci-wave gap symmetry by setting Tij = +1(— 1) for j = i±x (j = i±y) 
respectively, and with x, y being the square-lattice basis vectors. The (i-wave solution (of 
the dx2-y2 form) is taken throughout in the following analysis, as it is the most favorable 



energetically (cf. e.g. Ref. l38l ) The superconducting order parameter can be rewritten in 



terms of the singlet and the staggered 7r-triplet components, namely 

Aij = (Cj4Cit)o = ^ ((Cj^Cit + CjtCii)o + {CjiC,^ - Cj^Cii)o^ 



with 



Alf = lr,,(A^ + A5), (10) 

Ajf = ir,,(A^-AB). (11) 

The superconducting order parameter Aij is defined on bond (ij) (nearest-neighbor pair of 
sites). To define the gap per site, we make use of the standard^ relation for d-wave solution 



= - ^jm) = li^A + Ab), (12) 

^^^^ = i E ^^^•(^^.^•(0 + = l(^A- ^b) e^^^ (13) 

where denotes the nearest neighbors of site i. The existence of the triplet component 
is inevitable even if there is no triplet channel in the pairing potential. Namely, the triplet 
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FIG. 1. (Color online) Spin-majority (blue, bigger arrows) and spin-minority (red, smaller arrows) 
electron spins in a system with the AF order and real-space superconducting gaps. binds 
two spin-majority electrons, and binds two spin-minority electrons and therefore, there is a 
priori no reason for these two gaps to coincide (as would be the case for no staggered vr-triplet 
component). In other words, the two distinct gaps make effectively the t ~ 4- i ~ t pairing 
components of the opposite-spin pairs distinguishable. 

component is dynamically induced by the singlet gap and antiferromagnetism.— From 
a microscopic point of view, this is also not surprising at all (see Fig. [1] for an intuitive 
illustration). An interesting feature of the superconducting gap defined by Eq. iQ is the 
nonzero momentum of Cooper pairs for the triplet component (it results from the e**"^""' 
term, in an analogy to center-of-mass momentum Q in the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) phase^"— ). A superconducting state with nonzero momentum has been investigated 
in a number of cases,— ^^^^^ even in zero external magnetic field.— The one presented 
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, except we consider both the microscopic t-J model 



here is analogous to that from Ref. 
and the limit of strong correlations. 

With the above assumptions, we can express the ground-state energy W [cf. Eq. (E])] 
as a function of the band filling n, the magnetization components and rriAF, the 

superconducting gaps and A^, and the hopping amplitudes Xija = (cLcjcr)o- We assume 
as nonzero first- and second-nearest-neighbor hopping integrals t and t', which yields 6 
different hopping amplitudes 



Xija = (cLcj<T)o e {xab^, Xabi, XAAt, XAAi, Xbbi^, Xbbi} ■ 



(14) 



The resulting expression for W is quite lengthy and has been presented in Appendix |X1 Next, 
ihe method proposed earlier—"— (our present formulation is analogous to that from 
481 1 . to solve the model in a statistically-consistent manner, we impose additionally 



as m t 
Ref. 



the constraints on all introduced mean fields by means of the Lagrange multipliers method. 



In effect, to carry out the subsequent analysis, we use the following energy operator 

ija ij 
-^^t\f^i<T -nta) + W{n,mAF,mFM,^A,^B,Xija) - fi^Ui^. (15) 



This method of approach is equivalent in the T — t- oo (/3 — > 0) limit to that presented 



in Refs. 
Eqs. 



52 and 
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Ref. 



- fl37l) with the corresponding equations from Refs. 



52 and 
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The equivalence can be seen from the comparison of Eq. f[T5|) and 

(e.g. Eq. (13) from 



521 provides an effective Hamiltonian with the operator part equivalent to our K). The 



Lagrange multipliers A-J^, Xif\ and Aj-^'' have the same symmetries, as the corresponding 
mean fields Xija, ^ij, and rijo-. We also assume they are spatially homogeneous. Namely, 



^iJ = A^i + cr Xmpni + cr XmAF ^ 



A^f + A!r)e^Q-, 



with 



AS) 



AT) 



Aas). 



(16) 
(17) 

(18) 
(19) 



After performing Fourier transformation of the operator part of K we obtain 
^ = Y1 ic^icMk^k + A(/i + A„ - A„,,J + W{A) 
+A JnA„ + mpM^mFM + rnAF^mAF + 4(AaAa^ 

+ 4y^(2XABaA^^g^ + XAA^XAAa + XBBa^BBa] 



AbAa, 



(20) 



where the primed summation runs over the folded (magnetic) Brillouin zone, by A we denote 
all the mean-fields, A is the total number of sites, and the four-component operator \E'j^ has 
the following components 



(21) 



The matrix Mk is given as 



-2A^^?7k Ck+Qt — 2A^^?7k+Q \ 



-^k4 -2Ai^)r/k 



kj. 



^kt 
Ckt 

V -2A^^?7k+Q C-k+Qi 2A^^?7k+Q -^_k+Qi / 



-2A^^?7k 'Ck+Qt 2A^^?7k+Q 



(22) 



where 



Cka = -O- (A 



'XAAa 




(23) 
(24) 
(25) 
(26) 
(27) 



^ko- — — — A„ — (tA, 



- 2ekA 



4(A 



+ A 



ek = 2 (cos /ca; + cos ky) , 
= 4 cos cos /cy, 



?7k = cos kx — cos /C; 



We have also used the fact that Xlk^k = Xlk^k ~ ^- Note that in the present formulation 
XmpM corresponds to sum of the magnetic field h and the correlation-induced field hco-r^^ 
(or equivalently the Lagrange multiplier /3 in the slave-boson theor j*^""— ) . Namely, Xmp^ = 
h + hcor = h + P, what is evident from comparison of Eq. flM|) with appropriate expressions 



Next, we determine the eigenvalues of Mk, as they correspond to quasiparticle excitations 
of the system. An analytic diagonalization of Mk produces very long expressions, and more 
importantly, expressions with square roots of possibly negative numbers. Therefore, having 
in mind their subsequent implementation to calculate the physical properties, we diagonalize 
this matrix numerically. Next, having determined the eigenvalues {-E'ki}i=i,2,3,4, we determine 
the generalized grand potential functional for the system of fermions, which is 



The physical (equilibrium) values of the mean fields and the Lagrange multipliers are ob- 
tained from the necessary conditions for J-" to have a minimum subject to the constraints, 
i.e.. 



taken from Refs. Mm 




^ = k,.=i,2,3,4ln(l + e-^^-) + A(/i + A„ - A„,,J + W{A) 



+A ^nXn + mpuXmpM + "^afA^^^ + 4(A^Aa^ + A^Aas) 

+4^(2XABaAx^a„ + XAAaX^AAa + XBBaXxBBa) ■ 



(28) 




(29) 



where by A we denote collectively the Lagrange multipliers. Equations dJ^ /dA = provide 
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the explicit analytic expressions for the Lagrange multipliers, i.e., 



A„ = -A-i(9„iy(l), (30) 

A^,,, = -K~'dm,,,W{A), (31) 

A^^, = -K~'d^^,W{A), (32) 

Aa, = -\K-^d^,W{A), (33) 

Aa, = -\^'^^^,W{A), (34) 

Kab. = -Ia^%,,W{A), (35) 

Kaa. = -\A''d^,,W{A), (36) 

Knn. = -\a-%,sW{A). (37) 

The above expressions can be utilized to eliminate Lagrange multipliers A from the solution 
procedure. Thus, we obtain 11 equations to be solved numerically for the mean fields A, 
instead of 22 equations for both A and A. The equations for the mean fields (obtained from 
dJ^/d\ = 0) have the following form 

= /3-'dxM\)-A{n-l), (38) 

= r'dx^,,,M\) - H^FM + 1), (39) 

= p-'d,^^J^iX)-AmAM, (40) 

= r'dx^J^iX) - AAAa, (41) 

= r'5A^,//3(A)-4AAs, (42) 

= r'dx^,,J^{\)-8AxAB., (43) 

= r'dx^^^J^iX) - AAxAAa, (44) 

= /3-'dx^,,J(si\) - ^AxBBa, (45) 



where 



M\)^ J2 Ml + e-/^^-)- (46) 



k,i=1..4 



The derivative dx^fp{\) is computed numerically with a 5-point stencil method (as it gives 
two-three orders of magnitude better precision than the standard 3-point stencil). For 
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example 



1 



12x 



fpi^n + 2x, XmFM, ^mAF, •••) + 8//3(An + X, XmFM, ^mAF, 
— 8//3(A„ — X, XmFM, XmAF, •••) + //3(A„, — 2x, XmFM, XmAF, ■■■ 

+0(x') 



(47) 



— * 

where we use the "equihbrium" values of A as given by Eqs. ( ]5U]) - (157I) . The step x is typically 
equal to x = 0.0001. Larger values of x would cause greater error in the above formula. 
Smaller values would result in a loss of numerical precision. We have verified that at /i = 
(where analytical formulas for the eigenvalues i?kj are available) the numerical computation 
of the derivatives according to Eq. fHTj) with the chosen step x = 0.0001 introduces error 
smaller than the precision of the procedure of solving the set of Eqs. (I38l) - (l45l) . The value of 
the step X has been chosen after an analysis of the error aX h = and the numerical-precision 
loss (for very small x the numerical-precision loss lead to impossibility of solving the set of 
Eqs. (I38l) - (l45ll with the given precision). 



III. RESULTS AND PHYSICAL DISCUSSION 

The equations ( l38l) -( l45|) are solved numerically with the use of GNU Scientific Library 
(GSL)^ on a grid of size A = 256 x 256. We use the gsl_multiroot_f solver_hybrids 
solver which implements the hybrids algorithm. We use the precision epsabs = 10"''. 
Namely, the procedure converges when the relation < epsabs is fulfilled (where 

the sum is taken over all equations, which have been brought to the form /« = and divided 
by A to ensure lattice- independent convergence conditions). We assume the following values 
of parameters: t = 3, t' = t/4 = 0.75, J = 1, and /3 = 500, what yields the temperature 
T = 1//3 = 0.002 ^ 0. In Table I the exemplary numerical values of the parameters have 
been provided for the sake of completeness. Numerical accuracy is on the level of the last 
digit specified. The energy scale has been set by taking the value of the exchange integral 



as unit, J = 1. For more details on the numerical procedure see Chapter 6 of Ref. 
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Table I. Equilibrium values of mean-field variables, Lagrange multipliers, 


free energy F and g 


rand potential functional J-" at /i = 0.3 and /3 = 500. 


V di Id UlC 


V dlU." 


V dl IdUlC 


\/ n In *3 
V dlLlC 


IJ 


3 2997360 




-5 1331996 




0000000 




3000001 


Tfl A TP 


8100315 


X 

^mAF 


2 5817963 




0922998 


Aa 


2730140 




0479298 


Aa 


4395630 




1218625 




4258402 




0.1218625 




0.4258402 




-0.0167505 




-0.1031297 


XAAi 


0.0275895 


^XAAi 


-0.0120561 


Xbb^ 


0.0275895 


AxSBt 


-0.0120561 


Xbbi 


-0.0167505 




-0.1031297 


F/A 


-1.0110048 




-4.2117488 



A number of stable phases emerge as solutions of the equations, depending on the physical 
condition (n, K). As we work with constant number of particles n, the stable phase is the 
one with the lowest free energy, defined by 



F = J^o + /i?^A, 



(48) 



where all the optimal values of mean fields and Lagrange multipliers (i.e. those being solution 
to Eqs. fl30p - (HS]) ) are inserted in the functional T and /i is the chemical potential. 

The exemplary phase diagram on the band filling n - magnetic field h plane is exhibited 
in Fig. [21 It can be seen that antiferromagnetic phase is the predominant one in the low-field 
regime and above n = 0.8. For n > 0.935 antiferromagnetism coexists with superconduc- 
tivit y, w hat amounts to a phase with nonzero three order parameters (similarly, as in e.g. 



Ref. |25| ). In the low-n part of the phase diagram (for n < 0.8) the saturated ferromagnetic 
(SFM) phase with mpu = n becomes the stable state. This phase is stable even in the 
— 7- limit. This is an interesting result, which adds to the discussion of ferromagnetism 
in the t-J—^ model. There is also a number of papers (see e.g. Refs. 



48 and 



63l ) analyzing 



the t-J model ([T]) with the Gutzwiller-type of approach with the parameters in a similar 
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Phase Diagram: T = 1/p = 0.002, t'=t/4, J=t/3=1 




0.98 0.96 0.94 0.92 0.9 0.88 0.86 0.84 0.82 0.8 
Band filling, n 

FIG. 2. (Color online) Phase diagram on the band filling - magnetic field plane. The phases 
are labeled as follows: AF+SC - phase with coexisting superconductivity and antiferromagnetism, 
AF - antiferromagnetic phase, PP - polarized paramagnetic phase, SFM - saturated ferromagnetic 
phase (with mpM = "■)• No stable pure superconducting solution has been found. For a further 
analysis we restrict ourselves to n = 0.97 as marked by the dashed vertical line. 

range (i.e. with n < 0.8 and similar values of tij and J). Some of those papers disre- 
gard completely the Zeeman-term influence, and this omission is justified when applying 
the model to high-Tc superconductors, where orbital effects dominate over the Pauli mag- 
netism. We have shown, that even at /i = the system may be completely spin-polarized 
and therefore, the inclusion of ferromagnetic polarization is important in treatment of 
the t-J model. Finally, our phase diagram can be compared (although this is not a direct 
comparison, as even at /i = we have mpM 7^ for the AF phase) to that obtained by 
other Gutzwiller approximation scheme,— in which the coexisting phase was stable up to 
the doping 5 = 1— 71 = 0.1, and for higher doping levels the pure superconducting state was 
stable. In our approach the antiferromagnetic phase is stable for such dopings instead. We 
comment on the strong antiferromagnetism in the following analysis. 

At band filling n = 0.97 the phase diagram (or the phase sequence as a function of field 
h) resemble those observed recently in the heavy-fermion compounds CeCo(Ini_2;Cd3;)3^ at 

12 



o 

CD 



-0.7 
-0.8 



>^ -0.9 

D) 

03 
C 

LU 



-1 
-1.1 

-1.2 



Free energy at n = 0.97 



AF+SC 



AF+SC 
AF 
SC 
PP 



AF 



PP 
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Magnetic field, h 



1.2 



FIG. 3. (Color online) Free energy per site as a function of magnetic field for the specified selection 
of phases. Types of order are marked explicitly at the bottom. The "SC" phase is the pure 
superconducting phase (i.e. with mAF = 0), which obviously has higher energy than other phases 
and hence does not appear in the phase diagram. The vertical dashed lines mark the phase 
boundaries: the AF+SC - AF line marks a continuous transition, whereas that for the AF - PP is 
discontinuous, as one can see by looking at the behavior of the slope dF/dh. 

doping X = 0.0075 and CeRhSi^ at pressure p ^ 17 kbar.— Namely, in low magnetic fields 
a phase with coexisting antiferromagnetic and superconducting orders (AF+SC) is stable, 
whereas for higher magnetic fields a continuous transition to a pure antiferromagnetic (AF) 
phase takes place, followed by a discontinuous transition to the polarized paramagnetic (PP) 
phase. The phases appearing at this band filling (n = 0.97) are analyzed thus in detail in 
the following. 

In Fig. [3] we have plotted the free energy curves for a choice of possible a priori phases. It 
can be seen (also from the following Figures) that the transition AF+SC — ?■ AF is continuous, 
whereas the transition AF — )■ PP is of the first order. Also, pure superconducting (SC) 
solution is unstable, and this holds for other band fillings as well. It can be concluded 
from Fig. [3] that antiferromagnetism is the "dominating" phenomenon, since the energy 
gain from developing antiferromagnetic order (which can be seen from closer look at the 
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FIG. 4. (Color online) Staggered magnetic moment (a) and ferromagnetic spin polarization (b) 
for the selected phases. Obviously, the staggered moment of the SC and PP phases is 0, and has 
not been plotted in (a). The magnetic moment value is insensitive to the projection (i.e., it is the 
same in both the correlated l^*) and the uncorrelated |^o) states). 

difference {Fpp — F^f)) is much higher than the gain from developing superconducting 
order {Fpp — Fsc)- Moreover, the energy gain from developing AF order within SC phase 
{Fsc — Faf+sc) is much higher than that from developing SC order within the AF phase 
{Faf ~ Faf+sc)- This observation can be compared to the results of the Variational Monte 
Carlo method,— in which the d-wave solution is only slightly higher in energy than the 



coexisting phase (see Fig. 2 of Ref. l66l ). 

In Fig. m we exhibit the magnetic moment per site of the system for different phases. 
Namely, we plot the staggered magnetization uiaf and the ferromagnetic magnetization 
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1^FM (spin-polarization). The staggered magnetization is close to the saturation value of 
rriAF = n = 0.97. Such overestimation of the staggered magnetization value over the Vari- 
ational Monte Carlo results^i^ is also present in the slave-boson approach.— This is not 
surprising, as the method we use^"— is similar in structure (the Lagrange multipliers are 
present in both methods) to the slave-bosons approach (for the discussion of the equivalence 



for the paramagnetic state see Ref. |51| ). The obtained ferromagnetic spin-polarization for 
the pure AF phase is equal to mpM = 1 — n at all magnetic fields. Also, it can be seen that 
development of the SC order within the AF phase alters by a small amount the staggered 
magnetization rriAp, which drops by approximately 1% (see Fig. H^). 

In Fig. [5] various superconducting gaps are shown. Namely we exhibit both the "uncor- 
related" gap for the wave function |\l/o); as well as the gap for the correlated wave function 
\^), the latter defined by Eq. (jlj) and labeled as Ac. In the pure SC phase the sublattice 
gaps are equal (A^ = A^), what amounts to the absence of the triplet component. Note 
that although the uncorrelated gaps (A^, A^) are larger in the pure SC phase than in the 
AF+SC phase, the correlated gaps (Ac'^\ Ac^^) are much larger in the AF+SC phase than 
in the pure SC phase. This very important conclusion means that the presence of anti- 
ferromagnetism supports superconductivity in the present situation. The opposite is not 
true as the staggered moment is slightly larger in the AF phase than in the AF+SC phase. 
Finally, the renormalized gaps are more than an order of magnitude smaller than their bare 
(uncorrelated) correspondants. 

The picture with large antiferromagnetic magnetization niAF (Fig. Hj) and small super- 
conducting gap (Fig. [5]) is consistent with the energy curves displayed in Fig. [31 To shift 
the energy balance towards the SC phase one could either decrease t', or increase J . By 
doing that within a wide parameter margin, the antiferromagnetic phase still remains a pre- 
dominant phase. Other possibility to weaken antiferromagnetism is to include additionally 
the intersite attraction {V^^-^^ Uifij) in the starting Hamiltonian. This has been shown to 
stabilize the d-wave superconducting state^^ (see Ref. [sG for the expression for the average 
value of this term within the extended Gutzwiller scheme we use). The strong antiferro- 
magnetism may represent an apparent feature of the Gutzwiller scheme used,— in which 
magnetization is not changed by the projection, as follows directly from Eq. ([6]). 

Finally, in Fig.[6]we display the quasiparticle energies (the Slater subbands) for the phases 
discussed above for n = 0.97. The crossing of one of the bands with the zero energy line at 
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FIG. 5. (Color online) Superconducting gaps versus magnetic field for AF+SC and SC phases, 
(a) Aa and Ab gaps obtained for the uncorrelated wave function |^o); (b) and (c) gaps for the 
state specified by the correlated wave function (labeled as Ac), as given by Eq. (b) shows the 
sublattice-specific A^ and A^ gaps, and (c) shows the singlet and triplet components of the gap. 
Note that the superconducting gaps (Ac) are enhanced in the AF+SC state (with respect to that in 
the pure SC state, which is however an unstable state). Also, the singlet and the staggered-triplet 
components are almost equal in the correlated state. 
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FIG. 6. (Color online) Quasiparticle energies (bands) for phases obtained at n = 0.97: (a) AF+SC 
phase {h = 0), (b) AF+SC phase with nonzero spin-polarization (h > 0), (c) AF phase, and (d) 
PP phase (for a different path in the Brillouin zone). The full Brillouin zone (d) is spanned by the 
vertices (itTr, ivr), whereas the folded (magnetic) Brillouin zone (a-c), by (ibvr, 0) and (0,ib7r). The 
characteristic points are the following: F = (0,0), S = (7r/2,7r/2), Y = (0,7r), and M = (vr, vr). 
Note that in (a-c) two of the energies E'ki (those with E'^i > 0) describe quasiparticle states, 
and the other two represent quasihole states. Also, in (d) only two energies are displayed, as the 
Brillouin zone is not the folded (magnetic) one, but the full Brillouin zone, because there is no 
antiferromagnetism in this case. The red (dark) curve in (d) describes a quasiparticle with a =1, 
and the green curve describes a quasihole with a The fully gapped electronic structure in 
(a)-(c) is caused by the magnetic (renormalized Slater) gap appearance in the AF+SC and the AF 
phases. Energy scale is in units of J. 



the S point of the Brillouin zone in Figs. [6]dc means that the quasiparticle excitations will be 
spontaneously created (are gapless), a circumstance leading to a nonzero spin-polarization 
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(cf. Fig. [6]), similarly as in the situation for the FFLO state.— Note also that for h = the 
AF+SC electronic structure is gapful for the d-wave superconducting phase, because of the 
presence of the Slater-type (magnetic) splitting. This is not true anymore for h > 0.8 (cf. 
Figs. Etc and Fig. S^) when a uniform ferromagnetic component appears. Also, the bands 
are sizably wider in the PP state. 

IV. SUMMARY 

In summary, we have carried out a detailed analysis of the coexistence of antiferro- 
magnetism and superconductivity within a microscopic t-J model, with the Zeeman term 
included. The strong correlations were accounted for by means of the extended Gutzwiller 
projection method. Also, the constraints assuring the statistical-consistency of the results 
have been included. We have obtained the phase diagram on the band filling-magnetic field 
plane, in which for the band fillings in the range n ~ 0.935 — 0.970 and with the increasing 
magnetic field, a series of phase transitions takes place. Namely, the system evolves from 
the coexisting (AF+SC) phase, through the antiferromagnetic (AF) phase, to the polarized 
paramagnetic (PP) phase. Also, the onset of superconducting order reduces the AF order 
parameter. By contrast, the superconducting gaps are enhanced by the presence of the 
AF order. In the AF+SC phase there are two superconducting gaps of an almost equal 
amplitude: the singlet and the staggered-triplet components. These features reflect in a 
qualitative manner the experimental findings in the CeCo(lni_^.Cd2;)3^'i^ and CeRhSi^ 
heavy fermion systems, although our model is too simplified to be quantitatively related 
to such complex heavy fermion systems. Additionally, both antiferromagnetism and super- 
conductivity originate from the same electrons and are driven by the same kinetic-exchange 
interaction. Note that the real-space pairing is the pairing without "boson glue" , i.e. without 
paramagnons. It is the mechanism of pairing arising entirely from interelectron correlations 
which is particularly effective when renormalized hopping and exchange interaction are of 
comparable magnitude. 

As mentioned earlier, it would be very interesting to perform similar analysis within 
the Periodic Anderson Model, as this might allow for a comparison with the experiments 
[work along this line is in progress in our group^]. Also, testing other Gutzwiller schemes 
seems crucial to verify if the strong antiferromagnetism and the tendency towards saturated 
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ferromagnetism are only the characteristic features of the utihzed renormahzation scheme, 
or represent a universal tendency of the projected t-J model. For that purpose, the inclusion 
of realistic, orbitally degenerate / level structure, not just pseudospin doublet of Ce^"*", 
would be desirable. 

One should also note that the present approach includes the effect of applied magnetic 
field only via Zeeman term (the Pauli limit). For discussion of high-temperature supercon- 
ductivity for h > the orbital effects should be incorporated. 
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Appendix A: Explicit expression for W 



We provide here the expression for W = {'Htj)o- This expression can be divided into 
parts coming from different terms of Hamiltonian with Wt = Ylija'^iji^la^j^) = 
JT.{^j)i{S!Sp + {SfS^ + SfS])), as follows W = Wt + Wj - AhmpM- The expressions for 
Wj and Wt are given by 
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Wj = 2JA(- 4:{xABiXABt + A^Ab) 



^m\^ + {ml,, - (-2 + n)2)2 _ 2mi^(m^j,, + (-2 + n)2) 
+ -((-m^i;' + mFM){mAF + mpu) 
-(4(A^(-1 + niAF - mFM){-l + niAF + "^fm) 
x(2 + m^F - nipM -n){2 + tuaf + mpM - n) 
+x\b^O- + "^AF - mFM){-l + rriAF + mpM) 
X (2 + niAF + "mpM - n){-2 + tuaf - mpM + n) 

+XABi{-^ + ^AF - mFM){^ + rUAF + "^Fm) 

x(2 + m^F - - n){-2 + tuaf + "^fm + 

+A|(1 + rriAF - mFM){l + rriAF + mpM) 
x(-2 + niAF - mpM + n){-2 + tuaf + "T-fm + n))) 
/((2 + m^F - rnpM - n){2 + tuaf + - n) 

x(-2 + niAF - m-FM + 'n){-2 + tuaf + ^fm + "-)))) (Al) 



and 



... MABi^A^B^XAB^{'^x\Bi + T^AF-{'^ + mFM-n 



{m\p - (2 + mpM - nf ) \/-m\p + (-2 + m^^M + '^)^ 



-)t 



{XBB^{-4:xlBi + (-2 + rriAF - rriFM + nf)) 



\/~'^\f + (2 + vfiFM - nY{m\p - (-2 + mFM + "^Y 
+ (2(-l + n)( 



(2 + vriAF - mFM - n) 
XBBii-'^ + T^AF - mFM + n){-Ax\Bt + (^2 - mAF + mFM + 



-2 - niAF + "^FM + n 
I (-2 + m^i. - mpM + '^)^ 

+ (2(-l + n){xAA-\[-^X^AAi + (2 + m^i. + mpM - ^)^)(-2 + m^ir + vhfm + 

-Xaa4(2 + mAF + 'TiiT'M - '^)(-4xiAt + (~2 + "^af + mpM + nf))t') 

I ((2 + m^F + '^FM - '^)^(-2 + mAF + mpM + '^)^))- (A2) 
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